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Abstract 



We present a new general construction of recursion operator from zero curvature rep- 
resentation. Using it, we find a recursion operator for the stationary Nizhnik- Veselov- 
Novikov equation and present a few low order symmetries generated with the help of this 
£H . operator. 



In the present letter we suggest a new method of construction of recursion operator using 
the zero curvature representation. Unlike the majority of the hitherto known methods, see, 
e.g., |[|H|J and references therein, ours is immediately applicabile to both evolutionary and 



non-evolutionary systems, and gives not only the recursion operator, but also its inverse, lead- 
ing to the 'negative' part of the hierarchy of the system in question. We apply the method to 
the stationary Nizhnik- Veselov-Novikov (NVN) equation for which no recursion operator has 
■ been found so far. 

Let F = be a system of PDEs in two independent variables x, y for the unknown n-com- 
ponent vector function u = (u 1 , . . . ,u n ) T , where the superscript £ T" denotes matrix transposi- 
tion. Let this system have a zero curvature representation D y (A) — D X (B) + [A, B] = 0, where 
A and B take values in a (matrix) Lie algebra g and depend on A, x, y and u and its derivatives. 
Here D x , D y are the operators of total x- and y-derivatives, see, e.g., Ch 2 of 0, and ||. Note 
^ ■ that A and B may involve an essential (spectral) parameter A. 

Consider a (possibly vector or matrix) function P of x, y, u and its derivatives. Then 
the directional derivative of P along an n-component vector U = (U 1 , . . . , U n ) T is given by 

P'P] = El=iEZ=o( dP / 9u %M D i( Ua )> where u oo = u a , u% = ff+W/dxidyi, cf, e.g., 0. 
In P'[U] is called a linearization and denoted EpU. 

Let U be a symmetry of the system F — 0, that is, let U satisfy F'[U] = on the solution 
manifold of F = Consider a g- valued solution S of the system 

D X (S) — [A,S]=A = A'[U], D y (S) -[B,S]=B = B'[U]. (1) 

Assume that we have found n linear combinations U a = ^ • a a,!J Sij of entries Sij of S, 
a — 1, ... ,n, with the property that U = (U 1 , . . . , U n ) T is another symmetry of F — 0. Then 
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the linear operator D^o denned by U = 9{q(U) is a recursion operator, in Guthrie's || sense, 
for the equation F — 0. The coefficients a a,!J may depend on X,x,y and u and its derivatives. 

However, testing the above scheme on a number of known examples like KdV or Harry Dym 
equation shows that D\q generates the nonlocal symmetries that belong to the 'negative' part of 
the hierarchy of F = 0. Then we should, if possible, invert D\o in order to obtain a 'conventional' 
recursion operator £H, which will generate the 'positive', local part of the hierarchy in question. 
The inversion is an algorithmic process described in || . Note Q that if the coefficients of the 
recursion operator are local, then so are the coefficients of its inverse. 

Let us now apply this procedure to the stationary NVN equation 
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recently studied by Ferapontov |H| , see also Rogers and Schief |TT| , in connection with isother- 
mally asymptotic surfaces in projective differential geometry. 



The stationary NVN equation is a reduction of the NVN equation [^2| |13 
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obtained upon assuming that u,v,w are independent of t. The latter is well known to be 
integrable via the inverse scattering transform, as it has the Lax pair 



4>xy = U4>, lf>t = 



l/lyyy ~ 3V^ X + 3Wl/jy. 
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The NVN equation (|3]) is the first member of the hierarchy describing the deformations pre- 
serving the zero energy level of two-dimensional Schrodinger operator It also naturally 



arises in the theory of surfaces, see [|TTJ and references therein, and its modified version appears 
in the string theory 0, [15 . 

Upon setting |TT| ip = ipexp(Xt), where A is a constant, the Lax pair (|]) can be transformed 
into a zero-curvature representation for (Q) of the form D y (A) — D X (B) + [A, B] = 0. This 
representation involves an essential parameter A, and the matrices A and B belong to the 
semisimple Lie algebra sIq of traceless 6x6 matrices. They read 
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matrices with zero trace, and where A41 = —u yy + 3wu, B^ = —u xx + 3vu. 
Let U = (U, V, W) T be a symmetry of (@), i.e. let U, V, W satisfy 

D\U = D'lU + 3 [wDyU + u(D y W - D X V) - u x V + u y W + (w y - v x )U - vD x U] , 
D„U, D y V - n TT 



D X W 



D X U. 



(5) 
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Consider a traceless 6x6 matrix S that solves ([!]), where A, B are given by (|5|) and 
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are the directional derivatives of the matrices A and i? (|5D along the vector (U, V, W) T , An = 
-D 2 y (U) + 3wU + 3uW, £41 = -D 2 X {U) + 3vU + 3uV. 

The next step is to find linear combinations U,V,W of entries of 5* that solve (||). A 
straightforward but tiresome computation shows that for A and B (f|) these are U = — S35 + S 2 q, 
V = S54 + S12, = — Sis — 564, i.e. if U = (U, V, W) T is a symmetry of (0), then so is 
C7 = (£/, V", iy) T . Hence, the linear operator 9io mapping U to U is a recursion operator for (Q). 

However, the application of 9l to the simplest symmetries of (0), e.g. to the zero symmetry, 
yields nonlocal symmetries of (0), so we should invert 9^o in order to obtain a recursion operator 
91 = 1 generating hierarchies of local symmetries for (Q). 

It turns out that our 9^o is invertible only for A 7^ 0. Inverting 9to involves solving a system 
of algebraic and differential equations for the components of 91, which is a fairly tiresome but 
algorithmic process. For the sake of simplicity we set all the integration constants to zero. 
Then 91 = A9^ — |A 2 id, where id is the identity operator, is independent of A and provides a 
conventional recursion operator for (0). 

The action of 91 on a symmetry U 
Here Z = (Zi, Z 2 , Z 3 , Z 4 , Z 5 ) T is a general solution of the system 

D x Zi = U, D y Z 1 = W, D X Z 2 = V, D y Z 2 = U, 
D X Z 3 = D 2 y U - 3(Wu + Uw), D y Z 3 = D 2 X U - 3{Vu + Uv), 
D X Z 4 = -\vD 2 y U + \u x D y U + {-\u xy + u 2 + vw)U + {-\u yy + \ 

+2uw — v 2 )V + uvW + uu x Zi + {wu x + uu y — vv x )Z 2 , 
D V Z 4 = -\vD 2 x U + \v x D x U + {-\u yv + 2mu)C/ + uv V + 



(U, V, W) T of (D is given by 9t(l7) = £(17) + TIZ. 
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= -\wD 2 x U + \u y D x U + + u 2 + in«)Ef + mu;^ 

+ ( — | M xa: + 1^2/2/ + 2w — W 2 )!^ + («W X + t>U y — WW y )Zi + UU y Z 2 



Note that this system is compatible if and only if C7 solves 
The operators £ and 971 are of the form 
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£ n = D% - 6vD% - §uD 2 x D 2 - 15v x D s x - lu y D 2 x D y - §u x D x D 2 y + {-\ Uyy - 18v xx 
+ fuw + 9v 2 )D 2 + 9u 2 D x D y + {-\u xx + fuv)D 2 y + (-3^ - 12^ 
+ fwu x + 26uu y + 27vv x )D x + (26uu x + \vu y )D y - 3u xxyy - 3v xxxx 
+14wu xx + 20uu xy + hvuy y + 9vv xx + Y u xUy + 9v 2 — 4-u 3 — 28uvw, 

£ 12 = -f uD x - l -fu x D 3 x + (-fu xx + fuv)D 2 x + (-fu xxx + 7 -£vu x + 18uv x )D x 
-Quxxxx + I9vu xx + Auuyy + lOuv xx + ™u x v x + \u\ - Y2u 2 w - 4uv 2 , 

£l3 = -\uDy- + \UyD Z y + {-\Uyy + \uw) + (\U XXX ~ VU X - ^WUy ~ UV X + UWy)Dy 

+ 13UU XX + WUyy + UWyy + ~ \UyWy ~ \2U" 'v ~ 4:UVJ 2 , 

£ 21 = §DiD 2 y - fwDj - QuDlDy - fvD 2 x D 2 y - fu y D x - ^fu x D 2 x D y - fv x D x D 2 
H-W^v + f u " + fvw)D 2 x + {-^fu xx + fuv)D x D y + (-f u yy 
-fv xx + fuw + lv 2 )D 2 y + (-Uu xxy + fuu x + fvuy + fwv x )D x 
+ (-fuxxx + y VUx + f uv x)D y - 6u xxxy + 20uu xx + fvu xy + f wu yy 
+ yw% + fu 2 x + I9u y v x - 12u 2 v - 16uw 2 - 4v 2 w, 

£22 = ~±D X + \vD% + v x Dl + (-fuyy + fv xx + fuw - 3v 2 )D 2 x + {-fu xyy 

~^~~qV xxx ~\- ~^WU X -\~ -^UUy 7vV x ^jD x ~~g~Uxxyy g^xxxx ~^^U XX 

+Y uu xy + Y vu yy ~ f vv xx + Y UxU y ~ ^ v x ~ ^w 3 ~~ 16uvw + 4v 3 , 
£23 = -%UxDy + {\u xy - \u 2 )D 2 + (-fuxw + \wu x + 2uu y )D y - fu xxxx + y vu xx 
—wu xy + Y uu yy + y UVxx ~T u x v x + \u x w y — 12u 2 w — 4uv 2 , 

£31 = - SuD x - f vDlDy - lwD 2 x D 2 y - lfu x D\ + (-f v x - \w y )D 2 x D y 

-§UyD x D 2 + {-^fu xx + f uv + lw 2 )D 2 x + (-fuyy - fv xx + fuw 
+ fv 2 )D x D y + (-lu xy + |m 2 + \vw)D 2 y + (-±fu xxx + ±fvu x + f wu y 



+ fuV X + UWy)D x + (~U X yy ~ fv xxx + |wtl X + ^Uy + ^ 



;j ' ' .r 

~\~ TrVWy^ -Dy GUxXXX ~t~ ~3~VU XX ~\~ —WU X y ~\~ -^UUyy ~\~ ~T^UV XX ~\~ ~U X V X 

+ \u x w y + 7u 2 — \2u 2 w — 16uv 2 — 4vw 2 , 
£32 = -f u y Dl + (-9u xy + f u 2 )D 2 x + {-lQu xxy + fuu x + Uvu y )D x - fu xxxy 
+ Y uu xx + Y VUx y + i wu yy + Y u x + Y u y Vx + ^ u y w y ~ 12-u 2 f — 4-uw 2 , 

^33 = ~^f D y + f^-Dy + Wy D y + - 3U> 2 ) J D 2 . + (fWyyy ~ 7WWy)Dy ~ \U XX yy 

~~^~Q^yyyy ~^^^xx ~^UU X y -\~ AvUyy —WWy y -\- ~r^U x Uy 

-4-u 3 - IQuvw + 4w 3 ; 

^xxxyy ^1UU XXX ~\~ ~~^UU XX y ~\~ 3vU X yy -\- ^mJJyyy ~^"^yU xx ~\~ 9u x U X y -\- 3l) x Uyy 

+ WyUyy ~ ^UyWyy ~ 4:U 2 U X ~ 6VWU X ~ 8UVUy ~ 6UWV X ~ 4:UWWy, 

07^12 Hxxxxx ~t~ <-)VU XX x ~t~ ~^UU X yy ~\~ ~^UV XXX ~\~ \Qv x U xx ~t~ ^^x^lyy ~t~ ~^x^xx 8UWU X 

—6v 2 u x — 4u 2 Uy — 10uvv x , 
9^13 = ~\u xxx + 2vu x + 2uv x , 9Jl u = -2u x , = -2u y , 

07121 HxXXXy ^ > U'U'XXX 3VU XX y ~\~ 2WU X yy ~\- ~~^U X U XX QV X U X y WyU X y ~\- UyUyy 

+3u y v xx + \u x w yy — 9uvu x — flu)" 2 'u x — 9uwu y — 5u 2 v x + u 2 w y , 

27^22 g~'U XXX yy ~\- gV XXXXX ~\~ ~~^WU XXX ~\~ 5UU XX y ~\~ ~^VU X yy T^VV XXX ~\~ ~^U y U XX ~\~ 1 Oil X U X y 

+ fv x u yy — \v x v xx — 10vwu x — 4u 2 u x — lOuvUy + 4v 2 v x — 8uwv x , 
9^23 = - jU xyy + 2wu x + 2uu y , Wl 2 4 = ~2v x , Wl 2 5 = ~2ux, 

07^31 g~^xxxxx ~t~ Q^yyyyy ~%~VU XXX ~\~ -^UU X yy ~\~ ~^~UV XXX ^^^yyy ^-^^ X U XX ~\~ -^U x Uyy 

+ I0u x v xx — ^WyWyy — 10v 2 u x — 10uwu x — A^Uy — 1 Ouv v x + 4w 2 Wy + 2uvw y , 

O7T32 ^xxxxy ~t~ ^^^xxx ~t~ <-)VU XX y ~\~ \Ou x U xx ~\~ &V x U X y ~\~ "^^y^yy ~t~ ~UyV xx 

— \huvu x — 6v 2 u y — 3uwu y — 5u 2 v x + u 2 w y , 
9^33 = -\u xxy + 2uu x + 2vu y , 97l 34 = -2u y , 9Jl 35 = -2w y . 
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Note that the above formula £H(E7) = £(U) + defines a recursion operator in the sense 
of Guthrie ||, and the system (^) defines a covering over (§). Formally, we could express Zj 
from (0) as Zi = D~ l U, Z 2 = D~ l V etc., and thus write 91 as an integro-differential operator, 
as it became traditional in the literature, see, e.g., || |5|, [Tj. However, if we drop the y-part of 
(0), we encounter certain difficulties in constructing new symmetries, cf, e.g., [Tj], ffHj. 

As integrating (^) involves arbitrary constants, we have 91(0) = C\S\ + c 2 S 2 + ■ • • + C5S5, 
where q are constants, and Si, . . . , S 5 are symmetries of @ of the following form: 
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The repeated application of 91 to S\, . . . , S5 produces five hierarchies of symmetries of the 
stationary NVN equation (Q), which can be visualized as follows (numbers in the top line denote 
the orders of symmetries): 
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We conjecture that all these symmetries are local and commute, as it is the case for the sym- 
metries of orders 1, 3, 5, ... , 11. 

Note that (0) has a scaling symmetry S = xu x + yu y + 2u. The application of 9^ to S yields 
a nonlocal symmetry of seventh order, which we conjecture to be a master symmetry for (fj), 
meaning that commuting 9i(S) with any symmetry belonging to one of the five hierarchies, 
described above, yields (up to a constant multiplier) the next member of the same hierarchy. 
The repeated application of 91 to S yields an infinite hierarchy of nonlocal symmetries for (Q). 

We believe that 9i is hereditary in the sense of [Tj| , but we have not yet checked this because 
of the huge amount of computations involved. 

As a final remark, let us mention the nonstandard structure of nonlocal terms of 9^1 in (0): 
they involve the derivatives of components of the symmetry, what is quite unusual, cf, e.g., [19 
for another example of this kind and |2(J for a comprehensive list of known today integrable 
systems in (1+1) dimensions and their recursion operators. 
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